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Marine and hydrokinetic technologies, which convert kinetic energy from currents in open-channel 
flows to electricity, require inflow characteristics (e.g. mean velocity and turbulence intensity profiles) for 
their siting, design, and evaluation. The present study reviews mean velocity and turbulence intensity 
profiles reported in the literature for open-channel flows to gain a better understanding of the range of 
current magnitudes and longitudinal turbulence intensities that these technologies may be exposed to. 
We compare 47 measured vertical profiles of mean current velocity and longitudinal turbulence intensity 
(normalized by the shear velocity) that have been reported for medium-large rivers, a large canal, and 
laboratory flumes with classical models developed for turbulent flat plate boundary layer flows. The 
comparison suggests that a power law (with exponent, 1/a= 1/6) and a semi-theoretical exponential 
decay model can be used to provide first-order approximations of the mean velocity and turbulence 
intensity profiles in rivers suitable for current energy conversion. Over the design life of a current energy 
converter, these models can be applied to examine the effects of large spatiotemporal variations of river 
flow depth on inflow conditions acting over the energy capture area. Significant engineering implications 
on current energy converter structural loads, annual energy production, and cost of energy arise due to 
these spatiotemporal variations in the mean velocity, turbulence intensity, hydrodynamic force, and 
available power over the energy capture area. 

Published by Elsevier Ltd. 
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1. Introduction 


energy capture area (ECA) of the CEC machine. Fig. 1 illustrates 
typical profiles of current velocity and turbulence intensity in open 


The siting and design of river current energy conversion (CEC) channel flows, and demonstrates how the inflow characteristics 
technologies requires an assessment of the spatiotemporal varia- vary over the ECA of the CEC machine. The average hydrodynamic 
tion in the current velocity and turbulence acting on the proposed force and available power estimates over a representative period 

of record are calculated as 
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and 


= 1 = 1 
eas Cra 
P=5xadxpxwW=5 


xaxpx(14+3h +y) xT (2) 


where p is the fluid density, a is the ECA, u is the instantaneous 
horizontal current velocity component, and y=u’3/o? is the 
skewness coefficient (which is negligible). The instantaneous 
horizontal current velocity component is assumed to be the 
predominant velocity component, which is perpendicular to the 
ECA at all times and is defined as the sum of the mean velocity and 
velocity fluctuation 


u=uU+u' (3) 
The turbulence intensity is defined as 


lu =oy/il (4) 


hydrokinetic 
ail turbines 


water depth, D 
mean velocity 


z centerline 
energy capture 
area (ECA) 


turbulence intensity 


Fig. 1. Typical profiles of velocity and turbulence intensity in an open channel flow. 
The figure illustrates two possible hydrokinetic turbine configurations: a surface 
deployed vertical axis turbine (a) and a bottom deployed tower-mounted hori- 
zontal-axis turbine (b). 
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where o, = Vu'u’ is the standard deviation or root-mean-square 
(RMS) velocity. Detailed derivations of Eqs. (1) and (2) are 
provided in Appendix 1. These equations show that accurate 
assessment of the average hydrodynamic force and available 
power requires resolution of the turbulent fluctuations and turbu- 
lence intensity as well as the mean velocity. For example, a 
turbulence intensity equal to 20% at hub height of the CEC 
machine can increase the hydrodynamic force and power by 4% 
and 12.8%, respectively. The effects of turbulence should therefore 
be considered in the structural design and energy production 
calculations of a CEC machine. 

In medium to large rivers that are desirable sites for commer- 
cial scale CEC development, defined here as open-channel flows 
with 50th percentile depths that exceed one meter and 50th 
percentile currents greater than 1 m/s, collecting accurate mea- 
surements of the instantaneous velocity u needed to calculate 
basic inflow metrics (such as T, Iu, F, and P) is challenging [1]. River 
flows also exhibit great temporal variability of discharge and depth 
over time scales varying from minutes to days. For example, the 
USGS gage data for daily discharge and stage on the Missouri River, 
illustrated in Fig. 2, shows that discharges and water stages can 
increase tenfold during extreme episodic flood events over an 
approximately thirty-year period. This large variability in river 
discharge and depth is commonly observed and will affect the 
magnitude and distribution of mean velocity and turbulence over 
the ECA of a CEC machine over its design life. However, profile 
measurements needed to obtain meaningful statistics on the 
spatiotemporal variability of such important inflow metrics would 
rarely if ever be available due to the high cost and difficulty in 
obtaining such measurements. To obtain such measurements, an 
acoustic Doppler current profiler (ADCP) is the most practical and 


si 


Jul98 


= 
o 


Gage height [m] 
N wo A a oO N œ o 


Jan04 Jul09 


4.8 


Gage height [m] 
a 


46 


44 


4.2 


4 1 1 
08/03/02 08/13/02 08/23/02 


09/02/02 
Date 


Fig. 2. Daily flow and gage height time-series record for approximately thirty-year period of record (POR) on the Missouri River, Nebraska (USGS 06610000). The inset plots 
show the flow and gage height time series during field measurements by Holmes and Garcia (2009). 
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economical instrument to measure the vertical profile of u at a 
given site. But unlike an acoustic Doppler velocimeter (ADV), an 
ADCP typically cannot resolve the higher frequency turbulence in 
the inertial subrange and exhibits greater instrument noise [2-4]; 
therefore, ADCPs usually cannot provide an accurate estimate of 
the turbulence intensity Iu. 

Alternatively, mean velocity and turbulence profiles in rivers may 
be simulated using numerical or analytical models. Analytical models 
include simple power and logarithmic laws for the vertical profiles of 
mean velocity for a flat plate turbulent boundary layer flow [5] and 
exponential decay models developed by Nezu and Nakagawa [6] for 
profiles of the longitudinal (predominant), transverse, and vertical 


root-mean-square velocities, wu’, Vvv, and /w'w’, respectively. 
An evaluation of flat plate boundary layer models by Nezu and 
Nakagawa [6], who compared model predictions of turbulence 
profiles with measurements from large rivers and canals reported 
by McQuivey [7], suggests that these models can perform reasonably 
well despite simplifying assumptions of steady uniform flow with 
fully-developed boundary layers. These assumptions are generally 
violated for most river flows where the channel geometry, rough- 
ness, mean-section depth, and bulk velocity typically vary along the 
longitudinal direction. River currents are often three-dimensional as 
a result of variations in river alignment and vortices forming at 
different length scales from bed forms and in-stream structures [8], e. 
g. bridge piers. Pressure gradients associated with nonuniform sur- 
face profiles and bed forms, including ripples, dunes and anti-dunes, 
and wind shear can cause significant departures from semi- 
theoretical models that estimate mean and root-mean-square velo- 
city profiles [9,10]. Nevertheless, the simple analytical models can be 
useful for obtaining first order approximations of inflow character- 
istics at a CEC site. 

In this study, our aim is to examine measured vertical profiles 
of velocity and turbulence in open channel flows to: 


(1) Assess the validity of classical models for predicting river 
inflow characteristics along a planned ECA of a CEC machine. 

(2) Determine the range of velocities and longitudinal turbulence 
intensities that CEC technologies may be exposed to over their 
design life as a result of the large variation of water depth 
typical for river flows. 


2. Methods 


While ADCP measurements in medium and large rivers are 
now reported with increasing regularity in the literature, data is 


usually limited to mean velocity measurements only, e.g. [11], [12]. 
This present study examines previously reported vertical profiles 
of predominant RMS velocity measurements (i.e. turbulence) as 
well as mean velocity. Published profiles from five rivers 
[7,13,15,16], a large water supply canal [7], and a laboratory flume 
[7] were examined. Aside from the flume measurements we 
limited our literature review to open channel flows which meet 
our criteria for “medium to large rivers” as described in Section 1 
since these are the sizes of channels which provide the most 
interest to CEC developers. The ranges of bulk flow properties 
reported during data collection, including discharge Q, water 
profile depth D, and section width W, are summarized in 
Table 1. Three of the rivers had depths of at least two meters with 
velocities exceeding 1 m/s at the time of measurement. For several 
sites, mean annual discharge Q,, records were available to com- 
pare with discharges reported during the time of measurements. 
All Reynolds numbers R were above 4 x 10°. Froude numbers F for 
all cases indicate subcritical flows with the largest F observed for 
the Hurunui River [15]. 

In order to derive a shear velocity from the reported field 
measurements, the longitudinal velocity profiles were compared 
with the log-law [6], which describes the vertical distribution of 
longitudinal velocity in open channel flows with rough bound- 
aries, 


z = (ie) 48.5 (5) 

We computed the shear velocity Us = ./7./p from the slope of 
the best fit line through the time averaged velocity profile. With 
the exception of the Holmes & Garcia [13] data set, many of the 
vertical profiles contained only one or two data points in the wall 
region (z/D < 0.2) where the log law applies. For these cases, data 
from the entire vertical profile was used in this fitting method to 
calculate the shear velocity. Results for shear velocity using this 
method compared well with the values reported by the original 
investigators. Profiles with significant non-monotonic behavior, 
due to 3D flow effects or vortex shedding, were eliminated. This 
non-monotonic behavior was identified when the fitting proce- 
dure for the shear velocity resulted in a poor coefficient of 
determination, R?. Any profiles with R? values less than 0.8 were 
discarded. Of the velocity profiles reviewed, only 12% were 
eliminated because they exhibited non-monotonic behavior or 
could not be well represented by the log law. The final dataset 
included 47 profiles. 

Measured profiles of u non-dimensionalized Umax by were 
plotted to compare with the power law equation used for fully 


Table 1 

Bulk flow properties of reviewed open channel flow data. 
Investigators Site Qn Qt pt we Re® Fre 

(m?/s) (m?/s) (m) (m) (10°) 

McQuivey [6] Mississippi 19000°* 7900-9200 74-16 570-890 3-9 0.06-0.17 
McQuivey [6] Missouri 910° 890-920 2.9-3.1 200-210 4-38 0.19-0.35 
Holmes & Garcia [12] Missouri 2200° 1400 49 350-400 5-9 0.13-0.17 
McQuivey [6] Rio Grande canal NR 14-26 0.85-0.91 21-22 0.8-1.3 0.36-0.49 
Nikora & Smart [14] Hurunui NR 250 1.1-1.2 85-90 1-5 0.70-0.79 
Carling et al. [10] Severn NR 100 NR NR 3-6 0.10-0.16 
McQuivey [6] 2.44 m flume NR 1-2 0.33-0.53 2.44 0.4-0.8 0.69-0.74 
Gunawan et al. [15] Blackwater NR 1.87-3.57 0.65-0.97 5.74-6.75 0.2-0.3 0.14-0.19 


NR=not reported. 


a Mean annual discharge from nearest USGS station #07289000, record period: 2009. 
> Mean annual discharge from nearest USGS station #06610000, record period: 1953-2009. 
€ Mean annual discharge from nearest USGS station #06935965, record period: 2001-2010. 


4 Reported by original investigators at time of measurement. 
€ Derived by authors, using depth averaged velocity and local water depth. 
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developed flat plate boundary layer flows 


L-6" 9 


Umax 


and to determine the fitted power law exponent 1/a. The deriva- 
tion in Appendix 2, shows that a value of a=6 is consistent with 
the Manning's equation. 
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Fig. 3. Histogram of water depths at a planned turbine deployment site at the 
Mississippi river near Baton Rouge, LA. 
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Similarly, measured profiles of RMS velocities normalized by 
the shear velocity U= \/to/p were plotted to evaluate the accu- 
racy of the exponential decay models developed by Nezu and 
Nakagawa [6] for steady uniform flow in smooth laboratory flumes 


vV WU /Ux = 2.30exp(—z/D) (7a) 
V VV /Ux = 1.63exp(—z/D) (7b) 
Vw'w’ /Ux = 1.27exp(—z/D) (7c) 


These expressions, Eq. (7a)-(7c), are universal for smooth 
boundaries between (0.1 and 0.2) <z/D<0.9, independent of 


Reynolds and Froude numbers, and show that Vu'u'/u,> 


Vv'V'/Ux > \/ww'/u,. They do not apply near the wall, approxi- 
mately z/D < (0.1 to 0.2), as the no-slip condition requires turbu- 
lence intensities to decrease from a maximum value to zero at 
z/D=0. Nor do these equations apply in the free surface region 


above z/D < 0.9, where \/w'w’/u, is dampened. 

To demonstrate how to evaluate the spatiotemporal variability 
of mean velocity and turbulence acting on the ECA of a CEC 
machine, we performed a case study based upon conditions at 
Scotlandville Bend, Mississippi River near Baton Rouge, LA. This 
location is the site of a planned commercial river energy project by 
Free Flow Power (FFP) LLC, and was the basis for the river current 
hydrokinetic resource reference model proposed by Sandia 
National Laboratories [17]. This hydrokinetic resource reference 
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v R.G. canal: McQuivey (1973) dunes/flat 0.8-1.3 0.36-0.49 
> Hurunui: Nikora & Smart (1997) gravel 1.0-5.4 0.70-0.79 
Severn: Carling et al. (2002) gravel 0.9-6.4 0.10-0.22 
© 2.44m flume: McQuivey (1973) alluvial flat 0.4-0.8 0.69-0.74 
+ Blackwater: Gunawan et al. (2010) sand 0.6-1.0 0.24-0.29 


Fig. 4. (a) Mean longitudinal velocity profiles. (b) Longitudinal RMS velocity profiles. The dashed horizontal line indicates z=0.5 m. CEC technologies will operate at depths 


greater than 0.5 m off the bed. 
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model included an idealized trapezoidal river cross-section based 
on river transects with the low water reference plane (LWRP) at 
1.2 m NAVD (North American Vertical Datum), channel bottom at 
-11.5 m NAVD, and top of bank at 18.5 m NAVD. The FFP turbines 
are ducted axial-flow turbines that are three meters in diameter 
(d), which can be deployed from towers mounted to the river 
bottom or pontoons at the water surface. Although the US Coast 
Guard requires that the turbines be deployed at this site 6.1 m 
below the LWRP so as not to pose a hazard to river navigation, we 
performed our analysis neglecting this constraint—making this 
work applicable to other CEC technologies at possibly different 
sites where deployment depth restrictions may not apply. 

The variation of the mean velocity and turbulence acting on the 
ECA is examined at low and high water depths for two deployment 
scenarios: a turbine deployed from a tower with the rotor hub 1.5d 
above the river bed, and a turbine deployed from a pontoon with 
the rotor hub 1.5d below the water surface. Based on water stage 
measurements over a 25-year period of record from the nearest U.S. 
Army Corps of Engineers’ gage, located 8 km downstream of the 
site, the low water depth is defined for the purpose of this study as 
the 25th percentile value (D=14.74 m) and high water depth is 
defined as the 90th percentile value (D=21.54 m). The LWRP, at 
12.7 m, is the 2nd percentile water depth. The results of the 
statistical analysis of water depth are summarized in the frequency 
and cumulative frequency histograms shown in Fig. 3. 


3. Results and discussion 


All mean longitudinal velocity u profiles and RMS velocities are 
plotted in Fig. 4a and b, with the elevation above the bed z plotted 
on the logarithmic scale to accentuate the spread of depth and 
velocity in the data examined. We assume that z=0.5 m above the 
bed is a conservative lower limit on the location of the lower 
boundary of the ECA. Positioning the ECA closer to the bed would 
cause undesirable effects, including reducing the available power 
that can be harvested in the water column, increasing hydrody- 
namic load asymmetries and turbulence levels, and increasing the 
risk of strike or fouling by sediments or mobile bedforms. The 
mean currents u are observed to equal or exceed 0.5 m/s, with 
maximum velocities near the water surface and generally ranging 
from 1 to 3.7 m/s. Given that flow measurements for the Mis- 
sissippi River by McQuivey [7] were taken when the flow was well 
below the mean annual discharge Q,, (Table 1), one would expect 
higher maximum @ at higher z for flows with Q>Q,,. The 
corresponding longitudinal RMS velocity profiles, shown in 
Fig. 4b, increase exponentially from the free water surface to the 
near wall region. For fixed beds, the no-slip condition requires that 
the turbulence intensity and all components of the Reynolds stress 
tensor are zero at the bed, but rivers typically have mobile beds 
with a non-zero mean velocity and Reynolds stresses. The RMS 


velocities Vu'u’ are observed to range between 0.04 m/s and 
0.5 m/s. At high currents, such as the River Hurunui, which has a 


maximum velocity of nearly 4 m/s, the RMS velocityV/u'u’ is at a 


maximum of 0.5 m/s. Values of vwu for the Holmes and Garcia 
[13] data are noteworthy because they are quite high with respect 
to their mean velocity values. 

Longitudinal turbulence intensity I, profiles with z normalized 
by the flow depth D are shown in Fig. 5, and the root-mean-square 


velocity v u'u’ is plotted on a logarithmic scale to show the large 
variation of turbulence. The plot shows that the maximum values 
of I, typically occur very close to the bed. The longitudinal 
turbulence intensity I„ increases exponentially from the free water 
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Fig. 5. Longitudinal turbulence intensity profiles with z normalized by the flow 
depth and the turbulence intensity normalized by the local mean velocity and 
multiplied by 100 to determine the percentage of turbulence intensity with respect 
to the local mean velocity. 


0.9 + 
0.8 t 
0.7 F 
0.6 + 
Q m = 1/5.4 
~ 0.5 5 
N 
0.4 + 
0.3 + m=1/3 
\ 
0.2 + X 
0.1 F 
$ 
0 0.1 0.2 0.3 0.4 0.5 
Tif ü maz 
Watershed: Reference: Bedform: Re(10°): Fr: 
a Mississippi: McQuivey (1973) dunes 3.9-38.6 0.06-0.17 
o Missouri: McQuivey (1973) dunes/flat 3.0-9.1 0.19-0.35 
ô Missouri: Holmes & Garcia (2008) dunes 4.8-5.6 0.14-0.17 
v R.G. canal: McQuivey (1973) dunes/flat 0.8-1.3 0.36-0.49 
> Hurunui: Nikora & Smart (1997) gravel 1.0-5.4 0.70-0.79 
Severn: Carling et al. (2002) gravel 0.9-6.4 0.10-0.22 
© 2.44m flume: McQuivey (1973) alluvial flat 0.4-0.8 0.69-0.74 
+ Blackwater: Gunawan et al. (2010) sand 0.6-1.0 0.24-0.29 


Fig. 6. Power law velocity profiles with z normalized by D and w@ normalized by 
Umax. The solid black line represents the best fit of the power law with exponent 1/a 
through the data, and the resulting best fit a=5.4 (R* =0.999). The grey lines (left to 
right) represent the power law with exponent 1/3, 1/6 and 1/12, respectively. 


surface z/D=1 into the near wall region z/D < 0.1, but the large 


scatter indicates that vV ww does not scale well with U. 

Field measurements of ŭ non-dimensionalized by Umax are 
compared with the power law Eq. (6) in Fig. 6. Based on the 
power law assumption, Umax occurs at the surface (z/D=1), but the 
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measured data shows that Umax can occur beneath the surface due 
to wind, wave and three-dimensional flow effects. The power law 
exponent 1/a was observed to vary from 1/3 to 1/12 between 
individual profiles, with a best fit value of 1/5.4 through all the 
data, which is in good agreement with the theoretical value 1/6 
that was shown to be consistent with the Manning's 
equation (Appendix 2). Variation in the exponent can be attributed 
to a number of causes, including measurement error, pressure 
gradients, roughness and three-dimensional flow effects. 


Field measurements of normal stresses, e.g. vV u'u’, normalized 
by shear velocity ux=4/to/p with exponential decay models 


developed by Nezu and Nakagawa [6] are compared in Fig. 7. 


The value of v ww /u, is expected to reach maximum near the bed, 
then gradually decreasing to zero at the wall (bed) [14]. A peak 


value of v wu’ /u.=3.8 was observed in the near-wall region by 
Holmes and Garcia [6]. This large departure from the exponential 
decay model may have been due to turbulence generated by the 
presence of dunes. The gradual decrease to zero, however, cannot 
be observed because of lack of measurements close to the bed. 
Obtaining such data is difficult in the field because of limitations 
inherent in the instrument used for measurements. Field mea- 
surements are in reasonable agreement with the exponential 
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Fig. 7. Exponential decay law profiles by Nezu and Nakagawa (1993) compared to field measurements, with z normalized by D and normal stresses, e.g. v uw normalized by 


shear velocity us = \/to/p. 
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decay models developed from laboratory flumes, although there is 
considerable scatter; likely due to a combination of measurement 
error and complex hydrodynamic effects summarized previously. 
The measurements by Holmes and Garcia [13] are the only 
measurements known by the authors to include normal Reynolds 


stresses ,/v'v'/u, and ,/w’w’/u, in large rivers (depths > 1 m and 
currents > 1 m/s). These turbulence measurements are in fair 
agreement with the exponential decay models, except near the 
surface where the models underestimate the data. The discrepan- 
cies can be attributed to measurement errors due to wave motion, 
wind shear effects, or close proximity of the measurement sample 
volume from the measurement boat. 

The effect of river depth (D) variation on the vertical location 
and size of the ECA relative to the characteristic profiles for the 
velocity and turbulence is shown in Fig. 8 for the geometry and 
flow properties described for the Scotlandville Bend site (Section 
2). The location and height of the ECA, like the mean velocity and 
turbulence profiles, are non-dimensionalized by the water depth 
D, which results in variation of the relative location and height of 
the ECA with water depth. This non-dimensionalization allows a 
comparison of the relative location and relative size of the ECA for 
extreme low and high water depth over the velocity and turbu- 
lence profiles for the near-bed and near-surface deployment 
scenarios. Furthermore, one can apply the power law (Eq. (6) with 
an exponent 1/a= 1/6) and the exponential decay model for the 
longitudinal turbulence intensity (Eq. (7a)) to quantify the effect of 
depth variation on the inflow metrics U, I,, F, and P between the 
low and high depth extremes—values for these inflow metrics 
based upon this methodology are given in Table 2. 

River depth variability will adjust the relative location and 
length of the ECA at the low and high water depth extremes as 
shown in Fig. 8 for the two deployment scenarios considered. 


non-dimensional water depth (z/D) 


1a — bed deployed — low depth 2a — surface deployed — low depth 


1b — bed deployed — high depth 2b — surface deployed — high depth 


Fig. 8. Effects of large depth variability on the location of the swept area (energy 
capture area) relative to the velocity and turbulence profiles. 


Table 2 


The high depth in this example is 46% greater than the low depth. 
The depth-averaged velocity U and maximum surface velocities 
Tmax (based on rating curves provided for the USACE gage site) are 
greater by 85% at the high depth compared to the low depth 
(Table 2). For the near-bed deployment scenario, the turbine hub 
height moves downward by 0.1(z/D) from the low and high water 
extremes. Application of the power law indicates that the mean 
velocity u at hub height (i.e. centerline of ECA) increases by 72%, 
the turbulence intensity I„ (Eq. (4)) at hub height drops by 38%, the 
hydrodynamic force increases by 193%, and the available power 
increases by 373%. Minimum and maximum values of these para- 
meters (at each end of the ECA), and percent differences between 
these values at the low and high depths, are also provided. For the 
near-surface deployment scenario, the turbine hub height moves 
upward by 0.1(z/D) from the low and high water extremes (Table 2). 
The mean velocity u at hub height increases by 89%, the turbulence 
intensity I, (Eq. (4)) at hub height drops by 29%, the hydrodynamic 
force increases by 241%, and the available power increases by 493%. 

Decreases in turbulence intensity can be deceptive because it is 
normalized by the local mean velocity. Gunawan et al. [18] showed 
that turbulence at a hydrokinetic energy site increases with mean 
velocities, while at the same time turbulence intensities decreases. 
In both near-surface and bottom deployment scenarios, observed 
reductions in turbulence intensity Iu =o,/U actually result in 
increases in the turbulence. For the near-surface deployment 
scenarios at hub-height, while I, drops by 29% between cases 2a 
and 2b, turbulence actually increases by 37% (Table 2). Similarly for 
the bottom deployment scenarios at hub-height, while I, drops 
from 0.24 to 0.15 between cases 1a and 1b, turbulence actually 
increases by 5% from 0.21 to 0.22. This observation implies that 
fatigue load increases with mean velocities while at the same time 
turbulence intensities decreases. 


4. Conclusions 


Accounting for the reduced accuracy of velocity and turbulence 
measurements in large rivers and the uncertainty in calculation of 
the shear velocity, the power law (Eq. (6)) and exponential decay 
model (Eq. (7a)-(7c)) reasonably approximated the average dis- 
tribution of the local velocity and turbulence measurements 
reviewed. The results suggest that these models can be used for 
first-order approximations of the spatiotemporal variations of 
river inflow parameters arising from large depth variations, which 
can then be used to predict hydrodynamic forces and available 
power that a CEC will experience over its design lifetime. The 
scatter does indicate that more accurate estimates of velocity and 
turbulence will require advanced numerical modeling supported 
by site measurements of mean velocity and turbulence profiles. 
More field measurements in rivers using acoustic Doppler current 
profilers (ADCP) and acoustic Doppler velocimeters (ADV) are 
needed to expand the limited data presently available for large 


First order approximation of extreme river depth variation on inflow characteristics relevant to CEC machines at river hydrokinetic resource reference site. 


Case D Zhuv/D U Umax u (m/s) 
(m) (-) (m/s) (m/s) — 
Min Max hub 
la 14.74 0.31 0.92 1.04 0.8 0.9 0.86 
1b 21.54 0.21 1.7 1.92 1.38 1.55 1.48 
A% 46 -32 85 85 71 73 71 
2a 14.74 0.69 0.92 1.04 0.96 1 0.98 
2b 21.54 0.79 1.7 1.92 1.82 1.87 1.85 
A% 46 14 85 85 91 87 90 


Ou (m/s) lu (=) F P 

— — (kN) (kW) 
Min Max hub Min Max hub 

0.21 0.21 0.21 0.26 0.23 0.24 29, 2.6 
0.22 0.23 0.22 0.16 0.15 0.15 7.9 12.3 
6 12 8 -35 -35 -36 188 371 
0.31 0.30 0.30 0.32 0.3 0.31 3.7 43 
0.40 0.41 0.41 0.22 0.22 0.22 12.6 25.5 
30 37 34 -31 =27 -29 238 500 


1a. Bottom deployed low depth, 1b. Bottom deployed high depth, 2a. Surface deployed low depth, 2b. Surface deployed high depth. 
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rivers with fast currents. With limited site data, including a gage 
station with a daily or hourly water stage record, a stage-discharge 
rating curve, and a stage-area rating curve, we demonstrate how 
to apply these simple models to derive first-order approximations 
of the spatiotemporal variation of velocity and turbulence over the 
ECA of a CEC at a given river hydrokinetic energy site if three- 
dimensional effects are minimal, the channel is sufficiently wide, 
and the exponent for the power velocity law is reasonably 
accurate. The depth-averaged velocity U can be approximated 
using the equation in Appendix 2, the maximum velocity for the 
power law calculated as Umax = (7/6)U (Appendix 2), and the shear 
velocity us = ./gDS. 

A case study was performed to demonstrate the effects of 
extreme depth variation on velocity and turbulence profiles for a 
three-meter diameter ducted axial flow turbine at a river hydro- 
kinetic resource reference site, using the power law (with its 
exponent set to be consistent with the Manning's equation) and an 
exponential decay model. The analysis showed pronounced 
changes to q, Iu, F, and P between the low and high water depth 
extremes, which highlights the importance of modeling these 
spatiotemporal effects at river hydrokinetic energy sites. Signifi- 
cant engineering implications on CEC machine structural loads, 
annual energy production, and cost of energy arise due to these 
variations in the mean velocity, turbulence intensity, hydrody- 
namic force, and available power over the ECA. 


Notation 


Q 


energy capture area (m°) 

denominator in power law exponent (-) 

turbine diameter (m) 

flow depth (m) 

Froude Number (-) 

time-averaged hydrodynamic force (N) 

skewness coefficient (—) 

turbulence intensity of the streamwise velocity compo- 

nent (-) 

von Karman constant (—) 

roughness height (m) 

1/a, power law exponent (-) 

time-averaged power (W) 

discharge (m™3/s) 

mean annual discharge (m?/s) 

Reynolds Number (-) 

fluid density (kg/m?) 

standard deviation or root-mean-square velocity (m/s) 

boundary shear stress (N/m?) 

depth-averaged velocity (m/s) 

instantaneous longitudinal, lateral and vertical velocity 

(m/s) 

pi local time-averaged longitudinal velocity (m/s) 

Tmax maximum time-averaged longitudinal velocity (m/s) 

u’, v’, w' instantaneous longitudinal, lateral and vertical velocity 
fluctuation (m/s) 

Ux shear velocity (m/s) 

v kinematic viscosity (m?/s) 

W river width (m) 

Zz height above channel bed (m) 


eA DNX NNYAR 


STETE 
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Appendix 1. Derivation of Reynolds averaged hydrodynamic 
force and power 


For hydrodynamic force calculation in Eq. (1), using Reynolds 
decomposition, u? can be expressed as follows: 
u? = (U + u)? 
u? = (T + uN) + u’) 


u? =U + 2uu’ + u’? 


uw =T + 20w +u? 


w= +u’? 


u = (1+ Ryu 

Therefore the power includes an additional term that is 
proportional to the turbulence intensity squared. 

For available power calculation in Eq. (2), using Reynolds 
decomposition, u? can be expressed as follows: 
u? = (T + u’)? 
w=U+uya+u’) 


w= (i + 20u’ +u)\u+u’) 


w= + Uu + 207u' + 20u? + Gu? +u2 
w= + 30u + 30u + u’? 
w= + 3u7u' + 3uu2+u3 


W =T + 30u2 +u3 
wWi/2 3 
= Uu? u3\_ 
w—=—(14+3-—~4+ T 
ip P 


= u2 u’3.¢3\ _ 
wW=(14+35+54|7 
rms 


w = (1 +3Ê + yR 


where skewness, y = u? /o?. 


Appendix 2. Derivation of Power Law Exponent consistent 
with the Manning's Equation 


Given the power-law vertical mean velocity distribution in a 
very wide open channel: 


T o 
—= aA — 
Ux ks 


where U is the local mean velocity; u, is the shear velocity; a is a 
constant; z is the distance above the channel bed; k, is the 
equivalent sand-grain roughness; and 1/a is the exponent. 


V.S. Neary et al. / Renewable and Sustainable Energy Reviews 26 (2013) 437-445 445 


Integrate the point velocity distribution to find the mean velocity: 


1 7? z 1/a auy D 1/a+1 
u=5/, aus (Z) dz= Ga + DD (2) ks 
v= a(R) 

(1/a+ 1) ks 


The maximum velocity occurs at z=D, and is given by 


7 D 1/a 
Umax = du k 
s 


Therefore, 


Kio 
mT O 


Comparing the expressions for U and Umax, we have 


Rewrite the expression for mean velocity, substituting 
Ux = \/gDS, to obtain 


_ a/g pi/a+1/251/2 
(1/a + 1)ki/? 


With a value a=6, this equation for U agrees with Manning's 
and Strickler's equations in the exponents on D, ks, and S. The 
Manning's equation is therefore compatible with a 1/6 power law 
velocity distribution, and the maximum velocity at the surface can 
be related to the depth-averaged velocity U by the expression, 
Umax = (7/6)U. 
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